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Abstract 

Operation parameters of magnetic quantum cellular automata are evaluated for the purposes of 
reliable logic operation. The dynamics of the nanomagnets is simulated via the Landau-Lifshitz- 
Gilbert equations with a stochastic magnetic field corresponding to thermal fluctuations. It is 
found that in the macrospin approximation the switching speed does not change under scaling of 
both size and distances between nanomagnets. Thermal fluctuations put a limitation on the size 
of nanomagnets, since the gate error rate becomes excessive for nanomagnets smaller than 200nm 
at room temperature. 
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I. INTRODUCTION 



The success of computing in the past 40 years was based on scahng the complementary 
metal-oxide-semiconductor (CMOS) transistors to the nanoscale size As it is anticipated 
that this scaling will approach limits defined by the quantum theory and thermodynamics ^ , 
the search is on for alternative logic technologies js, 4|, which would be able to supplement 
CMOS and have certain advantages compared to it. One promising technology among them 
is spintronics and nanomagnetics j^. 

Magnetic Quantum Cellular Automata (MQCA) have been proposed as one of the types 
of spintronic logic. MQCA are based on bistable nanomagnet elements that can perform 
basic logic operations by means of magnetostatic interactions. Nanomagnets are typically 
arranged in the shape of crosses - majority gates. A majority gate has three inputs and one 
output. The output's logic state is determined by the 'majority voting' of the logic states of 
the inputs. This gate is naturally suited for the magnetic dipole-dipole interaction that is 
the basis of MQCA. It also allows us to perform AND and OR logical functions by fixing one 



of the inputs, and (in combination with the NOT element) it can 
logical operation. Another type of spintronics - domain wall logic 



3e used to perform any 
a] can also be rendered 



in the form of majority gates A chain of nanomagnets carrying the logic variables was 
demonstrated by Cowburn and Welland [8]. Later, a majority gate based on these principles 
was proposed and experimentally implemented [9|]- 

To be a viable alternative to CMOS logic, MQCA must show that they can achieve a 
better (or at least similar) performance level at least in one of the benchmarks, such as size, 
speed, switching energy, bit stability and scalability. Some of these issues have been studied 



through simulations 



m 



111 ]. In this paper, our goal is to estimate how far can we push 



the limits of MQCA performance for all the benchmarks presented above. To this end, we 
will analyze a simplified model of MQCA that captures the basic physical principles that 
govern its behavior. We pay a special attention to the limitation stemming from the thermal 
fiuctuations of the magnetization. 

The paper is organized as follows. In Section [III we show how the bit stability of an 
MQCA element puts a lower bound on its size. In Section UTTl we introduce a simple model 
of the MQCA dynamics and use it to simulate the behavior of an MQCA majority gate 
and study the speed of a signal propagating along a chain of nanomagnets. In Section 
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we discuss the relationship between MQCA initiahzation and its stabihty. In Section |V] we 
simulate the effects of thermal fluctuations and study their impact on the error rate of the 
majority gate. Finally, in Section IVTl we summarize our results and present our conclusions. 



II. BIT STABILITY AND MINIMUM SIZE 

Our first step will be to study what type of constraints bit stability imposes on the size 
of MQCA. The basic element of MQCA is a nanomagnet that is used to store a single bit of 
information. Usually the nanomagnets are elongated along some direction which determines 
the easy axis of magnetization due to shape anisotropy. This bit is represented by the 
magnetization direction of this nanomagnet: "0" for the magnetization 'pointing up', i.e., 
along positive easy axis, and "1" for the magnetization 'pointing down', i.e., along negative 
easy axis. We thus need to require these two configurations to be stable and separated by 
an energy barrier to prevent bit-flip errors. Even though material properties such as the 
uniaxial anisotropy can be exploited to produce such a bistable system, shape anisotropies 
are more advantageous to produce such a result, and most proposals of MQCA are essentially 
based on this idea. 

Our mathematical model is based on the free energy of a nanomagnet with uniform 
magnetization M. It includes contributions from the shape anisotropy, material anisotropy, 
and the energy in the external magnetic field (see [l^ for a derivation). 

E = Ki{l- {in.e,,isf)V + ^fioM^Vm.Af.m - fioM,Vm.He,t, (1) 

where m = is the normalized magnetization (note that |m| = 1); Mg is the saturation 
magnetization of the material; V is the volume of the nanomagnet; /iq is the permeability 
of vacuum; Ki is the uniaxial anisotropy of the material and Gaxis is a unit vector in the 
direction of the easy axis; A/" is the demagnetizing tensor; and Hext is the external field. 
The demagnetizing tensor can be diagonalized by finding its principal axes, and its diagonal 
elements are positive and satisfy + Ny + Nz = 1. We will consider that our nanomagnet 
is a rectangular prism whose symmetry axes are aligned with the cartesian axes. We will 
also assume that the easy axis of the crystalline uniaxial anisotropy is aligned with the y 



axis. The explicit expression for these demagnetizing factors can be found in ISj. 



Let us consider the case of a vanishing external field. If a, b and c are the dimensions of the 



3 



nanomagnet in the x, y and z directions, we will assume that b > a > c, which corresponds 
to a rectangular prism elongated in the y direction. This choice of proportions translates 
into an inverse ordering of the demagnetizing factors (A^^ > A^^ > Ny). This makes the z 
direction the least energetically favorable. It is easy to see that the energy is minimal when 
the magnetization points in the y direction, either up or down. These are the two stable 
states that encode a bit of information. Then the energy barrier between these two minima 
is smaller when we consider the magnetization to be in the x — y plane. To compute this 
energy barrier, we just need to evaluate ([1]) in the x and y directions and subtract them. 
Then we have 



AE = E{m = e,) - E{m = By) = ^fioM^V 



- {Ny 



(2) 



From this equation we can extract a few useful facts: (z) the energy scale is given by ^fioM^V; 
(ii) the energy barrier, and hence the energy dissipation, scales down with the volume of 
the nanomagnet; {in) the geometrical anisotropy can be used to control the height of the 
barrier; [iv) the uniaxial crystal anisotropy can be seen as a correction to the geometrical 
anisotropy. 

The height of the energy barrier will determine the stability of the information stored 
in the nanomagnets, and hence its bit stability. The thermal fluctuations will cause the 
direction of the magnetization to vary and with a certain probability to turn over 90 degrees 
- the direction of the energy saddle point. After that the magnetization will flip to the other 
energy minimum. In a simple model the probability of the nanomagnet's magnetization 
flipping its direction due to thermal noise is given by pfnp = exp(— Ai^/fc^T), where ks 
is the Boltzmann constant. Since we are interested in MQCA as an alternative to CMOS- 
based logic, it is natural to require this error probability to be at least of the same order as 
that of CMOS transistors, which is of the order of 10^^^. This corresponds to the condition 
AE/kBT > 40. For room temperature we have ksT ^ 0.026eV^, and so we need AE ^ leV 
or larger. The energy barrier height gives an approximate estimate of the energy that will be 
dissipated every time we switch the magnetization direction of a nanomagnet. The exception 
would be slow adiabatic switching regime [l4| which we do not consider here. 

The lower bound on the height of the energy barrier, coupled with equation ([2]) allows us 
to extract a lower bound on the size of the nanomagnets. Since the energy barrier depends 
on the volume of the nanomagnet, any lower bound on it will translate into a lower bound on 



the volume. Assuming that the geometrical anisotropy is due to a 2:1 aspect ratio between 
the length and width of the prism, we can plot the values of thickness and width that are 
required to obtain a leV energy barrier. In Figure 1, we present this plot for three different 
materials: permalloy, CoFeB and Fe (with saturation magnetizations equal to 800k A/ m, 
1180kA/m and 1750kA/m, respectively.) For example, in permalloy, we can see that for a 




5 10 15 20 25 30 

Width [nm] 



FIG. 1: Thickness vs. width for a nanomagnet with an energy barrier of leV. The length is taken 
to be twice the width. 

thickness of Qnm, the nanomagnet needs to have a 15nm width and a 30nm length. Clearly, 
there is an advantage for higher values of the saturation magnetization, since we can achieve 
the same energy barrier height with a smaller volume (see Eq. ([2]).) 

It can be argued that the very high bit stability we are requiring (error rate ~ 10^^'') 
might be appropriate for a memory device, but may not need to be that high for a logic 
device. For MQCA, we only need the nanomagnets to maintain their state only during 
the time it takes to perform a certain computation. We might be able to reduce the size 
even further if we somewhat relax the bit stability requirements. However, given that the 
dependence of the error probability with the energy barrier is exponential, a small reduction 
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in size can have a huge impact on the bit stabihty. We can illustrate this point by repeating 
the plot in Figure 1 for permalloy, but for different values of the error probability (Figure 
[2J) We can see the rapid increase of the error probability even for a modest reduction on the 
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FIG. 2: Thickness vs. width for different values of the error probability (plot corresponds to 
permalloy, and a 2:1 aspect ratio.) 

size of the nanomagnet. This shows that the lower limit on the size of MQCA is a rather 
strong one if we want to preserve bit stability. 




III. DYNAMICS AND SPEED 

To estimate the speed of MQCA-based logic devices we will simulate their behavior using 
the Landau-Lifshitz-Gilbert (LLG) equations jisl . Since we are only interested in an 
order of magnitude estimate, we will skip the detailed micromagnetic simulations that are 



usually discussed in the literature ITj], and instead work with a very simple model of the 
MQCA. We will model each nanomagnet as a macrospin, but we will include the effects 
of geometrical and crystalline anisotropics in the computation of the effective field. This 



approximation is equivalent to assuming the magnetization is uniform over the whole volume 
of a nanqmagnet at any time, and neglecting magnetic moments higher than the dipole 
We expect this approximation to improve for decreasing nanomagnet 



moment 
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size, since the exchange interaction tends to force the magnetization to be uniform on a 
length scale of about lOnm. From our discussion in the previous section, we are interested 
in nanomagnet sizes of the order of tens of nanometers, so we are not that far from that 
regime. In any case, we are interested in an upper bound for the speed of MQCA-based 
logic, and a full simulation will most likely produce a slower device. 
The LLG equations [lol, [2^ for the macrospin model are 

^ - -TTS^-'-' X - (T^M... X (M<.) X Hg,,. (3) 

where M*^*^ is the magnetization of the i*^^ nanomagnet, Hejf is the effective field at the 
position of the i^'^ nanomagnet, 7 = g\e\/2meC = 2.21 x 10^mA~^ sec~^ is the Lande factor, 
and a is the Gilbert damping constant, which depends on the material and the environment 
of the nanomagnet and typically has values in the range 0.001 — 0.1. The effective field 
includes the contributions of any external field, the nanomagnet self-field and the field due 
to the dipole-dipole interaction with other nanomagnets. 

= H« - Af . M» + Y: Cfo-)M(^-). (4) 

j 

In this expression we are assuming that all nanomagnets have the same shape, and hence 
the demagnetizing tensor A/" is the same for all nanomagnets. This term can also include 
the effects of uniaxial crystalline anisotropy if we redefine the corresponding demagnetizing 
factor Ny ^ Ny — , where y is the easy axis of the crystalline anisotropy. The last 
term on the RHS of (HI) represents the dipole-dipole interaction between nanomagnets, and 
the matrices C(jj) are coupling constants determined by their size and relative positions. If 
(x*^*-*, y^'^\ z*^*-*) are the coordinates of the i*^ nanomagnet, we define the coordinate differences 
for a pair of nanomagnets as dS^^'' = x^*^ — x^^\ d^^^ = y^*^ — y^^\ d^^^^ = z^'^^ — z^^\ and the 



distance between nanomagnet centers as d^^^^ = y {dx''^Y + {dy''^y + {dx^^Y- The coupling 
constant matrices are given by 

34^)4^) 3(4^))2 - 34^)4^) 

34^')4^') 34^)4^) 3(4^'))2 - J 



(5) 



In our case, since all nanomagnets will be in the (x, y) plane, this expression simplifies since 
di^j) = 0. An important fact about the matrices C(jj) is that they are dimensionless, and 
hence invariant under scaling of both the sizes of nanomagnets and the distances between 
nanomagnets. We will see that this property is preserved by the LLG equations in our 
model. 

To simplify the simulation and analysis it is useful to normalize the LLG equations. This 
is accomplished using the following definitions 

(i) M« 
m = -irr- 

^^^^ " M, 

t' = t(7M,), (6) 

where now all the quantities on the LHS of ([6]) are dimensionless (note that [jMs] = sec~^.) 
With these rescalings and using vector identities and the obvious fact that ^^^^.M*^*) = 
0, we can rewrite the normalized LLG equations in an implicit form that simplifies the 
implementation of the simulation, 

= -m« X hi}f + amW x (7) 

These equations have the property that the value of the magnetization is constant, 
|m*^*^(t')| = 1, Vt', and this feature must be preserved in the discretized numerical model. 



To do this we employ the mid-point method 21[ with which this constraint is automatically 
satisfied. 

To estimate the speed with which MQCA switch, we simulated the behavior of the ma- 
jority gate. Let us first briefiy review its operation. The nanomagnets forming the gate are 
arranged as seen in Figure [31 We also include three nanomagnets with fixed magnetization 
that are used to simulate the inputs of the gate. The nanomagnets that form the gate are 
initially magnetized in the x direction, and then are left to evolve driven by the magnetic 
dipole interaction. The magnetization of each nanomagnet will tend to align itself with 
the field produced by the other nanomagnets at its position. The geometric anisotropy will 
force the magnetization to lay in the y direction, and the infiuence of other nanomagnets 
will decide if it ends pointing up or down. The fields of the three inputs will add at the 
position of the central magnet and decide its direction of magnetization, hence computing 
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FIG. 3: Majority gate: the thick arrows represent nanomagnets with fixed magnetization that 
simulate inputs to the gate. The remaining nanomagnet ahgn their magnetization in order to 
minimize the energy of the system from an initial magnetization in the x (horizontal) direction. 
The output of the gates can be extracted from the magnetization of the "output" nanomagnet on 
the right. 



the majority of the input signals. Finally, this signal can be read on the output magnet. 
Note that a signal that propagates horizontally is inverted every time it is received by the 
next nanomagnet (due to the antiferromagnetic coupling). This does not affect the function 
of the gate, though this feature must be tracked in order to correctly interpret the output 
of any MQCA-based gate. 

Again, in order to extract numerical estimates from the simulation, we specified the 
properties of the material [Mg and Ki) to be those of permalloy. The value of the Gilbert 
damping constant did not have a big effect on the simulation when confined to the typical 
range 0.001 — 0.01. We found that the typical gate time, measured as the time it took the 
output to reach 90% of its final magnetization, was about 700ps. An interesting feature of our 
model is that the normalized equations ([7]) are invariant under changes of scale, which means 
that the gate time is independent of size. Even though this is only true in this simplified 
model, and making less approximations will likely break this invariance, whatever effects this 
may have on the the gate time will likely be of higher order. This is in contrast to CMOS 
logic pj] as well as MQCA based on magnetic wires (rather than discrete nanomagnets) 

From the form of the normalized equations we can see that the speed of this gate will 



9 



depend on the material properties. In particular, the speed of the gate increases linearly 
with the saturation magnetization of the material. This follows from the scaling of actual 
time t with respect to normalized time t', as defined by 

Another issue that needs to be considered when analyzing the speed of MQCA-based 
information processing, is the speed of propagation of information. In MQCA this is ac- 
complished by chains of nanomagnets that are initially magnetized in the x direction, and 
evolve according to the dipole-dipole interaction propagating a signal, as can be seen in 
Figure H] for the case of a horizontal wire. Note that the antiferromagnetic coupling forces 

Propagating signal 



1 i 



FIG. 4: Signal propagation through a horizontal wire made up of a chain of nanomagnets. The 
antiferromagnetic coupling forces neighboring nanomagnets to become antiparallel. 



neighboring nanomagnets to be antiparallel. For vertical wires the coupling is ferromagnetic 
and the nanomagnets magnetization tends to become parallel. 

This evolution follows the same dynamical equations presented in the previous section, 
so we can use them to simulate the propagation of a signal along a chain of nanomagnets 
and estimate its speed. For nanomagnets made of permalloy with a width of about lOnm, 
separated by 15nm, the speed of signal propagation is around 100m/ sec, or equivalently, 
ISOpsec per magnet. This is of the order of the speed of sound, and would certainly limit 
the speed of an integrated MQCA chip if communication is done using the same principles 
as logic. This speed depends on the material through the saturation magnetization, but 
only linearly, so it is not likely that choosing a different material will solve this problem for 
MQCA. 
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IV. INITIALIZATION AND BIT STABILITY 



As discussed before, in order to run a MQCA-based logic gate it is necessary to initialize 
the magnetization of all nanomagnets in the x direction (i.e., the hard axis.) In terms 
of energy, this corresponds to placing all nanomagnets at the top of the energy barrier 
created by the geometrical and crystalline anisotropics (see Figure [5] a).) However, this 
configuration corresponds to an unstable equilibrium point for each nanomagnet, and it 
should be expected that small perturbations due to thermal effects and stray fields will 
randomly force the nanomagnets to relax to one of their stable configurations independent 
from the input signals. 

This is an important issue for any implementation of MQCA-based logic and some possible 
solutions have been suggested. One consists of exploiting the biaxial anisotropy of the 
material to create a stable configuration around the initialization direction, by generating a 



local minimum of the energy 



22|. If we consider the magnetization confined to the x — y 



plane, and note as 9 the angle between the magnetization direction and the x axis, the 
geometric and uniaxial anisotropy result in an energy profile proportional to cos^(6') as can 
be seen in Figure [S] a). The biaxial anisotropy introduces another term that is proportional 
to sin^(20), and by carefully choosing the parameters we can produce a local minimum for 
6^ = 0, as seen in Figure [5] b). 

This energy minimum provides a latch mechanism that keeps the initialized nanomagnets 
pointing in the x direction while the information from the input signal propagates through 
the chain of magnets. Once again, the effectiveness of this local minimum to trap the 
magnetization direction against thermal fluctuations, will depend on the height of the energy 
barrier around it (i.e., the energy difference between the peaks and the local minimum in 
Figure [5] b).) The reasoning of Section [TTl applies to estimate the energy of this barrier 
necessary to preserve the bit in its local energy minimum for sufficiently long time, and 
hence obtain an estimate of the strength of the required biaxial anisotropy. We realize that 
the requirements to the height of this barrier are contradictory - it should be high enough 
to prevent spontaneous transition to one of the global minima before the signal reaches the 
bit; it also needs to be low enough so that the signal can reliably switch it to the desired 
local minimum. In the next section we simulate the behavior of the majority gate, including 
the biaxial anisotropy, in the presence of thermal fluctuations. 
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FIG. 5: a) Energy profile for geometrical and uniaxial anisotropies. Stable configurations corre- 
spond to magnetization in the y direction (up or down). Magnetization in the x direction (initial 
configuration) is an unstable equilibrium point, b) Including a biaxial anisotropy produces local 
minima for magnetization in the x direction, stabilizing the initial configuration. 

V. THERMAL EFFECTS AND GATE ERROR PROBABILITY 



In this section we model the effects of the thermal fluctuations on the operation of MQCA. 
We especially focus on gate errors caused by spontaneous transitions from the local energy 
minimum after the initialization of elements of MQCA. 

Our simulations will use the stochastic LLG equations based on the midpoint rule derived 
by d'Aquino et al. in [21I. The only difference with the above model (Section IIII|) will be the 
inclusion of an extra term that represents the field generated by the biaxial anisotropy (we 
show in the appendix that the introduction of this term does not affect the useful properties 
of the discretized equations.) 

Let us start by considering the extra term in the normalized effective field that is respon- 
sible for the biaxial anisotropy acting on nanomagnet (i), 

2K, 



ef f (biaxial) 



(m«(l - (m«)2)x + m»(l - (m«)2)y + m«(l - (m»)2)i 
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The biaxial anisotropy constant K2 has dimensions of Jm~^. It is not difficult to show that, 
when restricted to the x — y plane, the contribution to the energy of this term is proportional 
to s\.v?{26), with 6 the angle between the magnetization direction and the x axis. In order 
to have a local minimum around = 0, the constant K2 must satisfy the condition 



K2 > K2min = \^0MlV 



(9) 



The thermal fluctuations manifest themselves as random variations of the overall mag- 
netization of the nanomagnet. We describe this process by the stochastic LLG equations 



23l . |2J|, which are obtained by adding a random force, or, in other words, a stochastic 
thermal magnetic field h!^\t) to the effective field in ([7]). Note that we are considering a 
different thermal field for each nanomagnet, since it is usually assumed that the thermal 
fluctuations in different nanomagnets are uncorrelated. The random thermal field h^\t) is 
assumed to be an isotropic vector Gaussian white- noise process with variance u"^, and so it 
can be expressed in terms of the Wiener process as h'^\t)dt = z/(iW^*\ Then, the stochastic 
LLG equations take the form 

(im« = -m« X (h^^ + h«^(„)) dt - m« x u dW« + a m« x cim«. (10) 

The value of u can be obtained from the fluctuation-dissipation theorem in thermal equilib- 



rium, and is given hj u = y^^y • 

Using ffTOj) . we simulated the behavior of the majority gate for various values of size, 
damping constant, and temperature. We fixed the saturation magnetization and uniaxial 
anisotropy to be those of permalloy, and studied the error rate of the gate as a function 
of K2 and for several values of the damping constant a. Starting with the nanomagnets 
initialized with magnetization in the x direction, each run simulated the evolution of the 
gate for 2000 ps. We considered the gate to be successful if the average of the output magnet 
during the last 300 was larger than 80% of the ideal output value (all runs used the same 
set of fixed inputs.) In any other case, we considered that the gate failed. For each value 
of the parameters K2 and a, we ran 1000 instances of the simulation. The results are 
presented in Figure O The error probability is plotted against the ratio of K2 to K2min, 
where K2min is the minimum value of the biaxial anisotropy that produces a local energy 
minimum around ^ = 0. If we increase we expect the error probability to decrease when 
we pass K2/K2min = I5 since the biaxial anisotropy becomes more effective in preventing a 
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FIG. 6: Error probability of the majority gate as a function of the (scaled) biaxial anisotropy for 
different values of the damping constant (T = 300K) . 

premature flipping of the nanomagnets spurred by the thermal fluctuations. On the other 
hand, if we increase the biaxial anisotropy too much, the local energy minimum is too deep 
for the signal to force the nanomagnet to flip. This is the behavior we can appreciate in 
Figure O For K2 > 2K2min, the gate becomes essentially frozen by the biaxial anisotropy; 
for K2min < K2 < 2K2min, the error probability seems to have a minimum for a certain 
value of K2, that depends on the damping constant. However, an important result of these 
simulations is that, for the particular temperature and size considered (SOram x 15nm x Qnm 
magnets), the gate error rate exceeds a certain minimum value, 15% in this case. The 
stabilizing effects of the biaxial anisotropy are either too weak, and spontaneous gate errors 
happen, or too strong, so that it prevents the normal evolution of the gate. 

One possible solution for the gate error probability will be to decrease the temperature. 
Then, thermal fluctuations will be weaker and smaller values of the biaxial anisotropy will 
be enough to keep the magnets magnetized in the x direction until the signal, in the form 
of the magnetization of a neighboring magnet in the y direction, reaches the magnet and 
makes it flip up or down. And since the required biaxial anisotropy is not too large, it does 
not freeze the magnet in its initial magnetization direction. We used our model to study the 
dependence of the gate error probability on the temperature, again running 1000 simulations 
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for each value of the temperature and the biaxial anisotropy, and then finding the minimum 
value of the error probability for each temperature. These results are presented in Figure [71 
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FIG. 7: Minimum error probability of the majority gate as a function of temperature. 

We can see that, as expected, the error probability decreases with decreasing temperature, 
although this decrease seems rather slow for temperatures above 150K. For temperatures 
below 307^, the error probability is below 0.001, but it could not be accurately estimated 
with the same number of simulation runs. 

Another approach to lowering the error probability of the gate is to increase the size 
of the magnets. We know that larger magnets have a larger energy barrier between the 
states of up and down magnetization. This increases the stability of the computational 
states of the magnets but it is not the reason why the majority gate becomes more reliable. 
The key parameter is the ratio of the height of the energy barrier surrounding the local 
energy minimum around the magnetization in the x direction, and the strength of the 
signal produced by neighboring magnets. We ran our simulations for different sizes of the 
nanomagnets, but keeping a 2:1 aspect ratio and a thickness of Qnm. Figure [8] shows these 
results. 

We can see that the error probability decreases fast with size. The mechanism for this 
behavior is the following. When we increase the size of the magnets following the pre- 
scription mentioned above, the depth of the local minimum increases, but this increase is 
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FIG. 8: Minimum error probability of the majority gate as a function of the length of the nano- 
magnets. 

approximately a linear function of the length. On the other hand, the volume of the magnet 
increase quadratically with the length (since we are keeping a fixed aspect ratio), and hence 
the strength of the magnetic field generated by the magnets also increases quadratically. In 
summary, the deeper local minimum does a better job stabilizing the magnet against thermal 
fluctuations, while the magnetic interaction grows faster, preventing the biaxial anisotropy 
from freezing the nanomagnets. From these results we see that nanomagnets with size less 
that 200nm have too high gate error probability and thus cannot be used to build MQCA. 

VI. SUMMARY AND DISCUSSION 

The goal of this work was to estimate the characteristics of an MQCA-based logic device, 
in particular the limits that can be achieved in terms of minimum size, gate switching time, 
switching energy, and gate error probability. To this end we analyzed a simplified model in an 
effort to understand how these features are affected by the basic parameters that characterize 
the MQCA. A reasonable requirement on the bit stability of these devices naturally leads 
to a lower bound on the size of the basic element of any MQCA. A nanomagnet must be 
at least 20nm long in one of its dimensions to prevent thermal fluctuations from inducing 
an error rate larger than that of today's CMOS transistors. Furthermore, reducing this size 
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results in a rapidly degrading bit stability of the components, making its applications in 
logic circuits less useful. Fault-tolerant design does not seem to help in this situation, since 
any reduction in the size of the nanomagnets will be offset by the increase in their number 
due to the overhead usually accompanies fault-tolerant implementations. Another way to 
push beyond this limit would be to work at much lower temperatures, but that regime will 
not be practical in the most common situations. 

The lower bound on size also provides us with an estimate of switching for MQCA. After 
initialization of an MQCA, energy is dissipated when the magnetization of each nanomagnet 
"rolls down" the energy barrier until it reaches a minimum energy configuration (like a ball 
rolling on curved surface in the presence of friction.) Then, the energy dissipated by each 
nanomagnet is just the energy it had at the top of the barrier, and that is just the height 
of the barrier. From the bit stability constraint we found that this height should be at least 
leV, and hence a MQCA could in principle dissipate about leV per nanomagnet. A logic 
gate such as the majority gate requires only five nanomagnets, so we could perform logic 
functions with a switching energy as low as a few electron-volts. This is a big advanta ge o f 
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MQCA over CMOS transistors, that requires several thousand electron-volt to operate 
This is, however, only a theoretical limit, and it does not take into account the practical 
difficulties of efficiently transferring such a small amount of energy to each nanomagnet. 

To estimate the speed of MQCA logic gates we considered a very simple model in which 
we approximated the nanomagnets by point dipoles when computing their interaction, but 
included the effects of geometrical and crystalline anisotropics through the computation of 
the effective field. This approach is less sophisticated than the micromagnetic simulations 
that have been used in the literature to study similar systems, but our goal was not to obtain 
a very detailed picture of the dynamics, but rather to have a good estimate of the fastest 
gate time MQCA can achieve. Our model includes all the fundamental elements of MQCA 
dynamics, and more refined simulations are likely to result in slower gate times. Using this 
simple model we found that the majority gate produces the required output in about 700 ps, 
which is slower than gate times expected from CMOS in the next few years. 

Another obstacle for implementing MQCA-based logic has to do with information trans- 
mission. In MQCA this is accomplished following the same basic principles as logic. Chains 
of nanomagnets propagate a signal through the dipole-dipole interaction. But the propaga- 
tion speed of this signal turns out to be around 100m/ sec, which is extremely slow when 
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compared with the speed of electric signals in a wire (typically around 10^ m/sec.) This is 
a huge disadvantage for any MQCA scheme. 

MQCA suffers from the problem that its nanomagnets are initialized in an unstable 
state before the computation. Thermal fluctuation will push the nanomagnets randomly 
into one of the stable states regardless of the value prescribed by the computation. It has 



been proposed [22] that exploiting the biaxial anisotropy of the material, can increase the 
robustness of the MQCA initial state against thermal fluctuations, preventing premature 
relaxation of the nanomagnets before the computation is complete. On the other hand, a 
strong biaxial anisotropy can completely freeze the dynamics, by trapping the magnetization 
in the local energy minimum of the initial state. We simulated the behavior of the majority 
gate in the presence of thermal fluctuations and analyzed the error rate of the majority gate 
for different values of the biaxial anisotropy, in order to find what are the optimal choices of 
the parameters. We found that for room temperature operation (T = 300K), the gate error 
rate has an impractically high value (> 1%) for all sizes of nanomagnet smaller than 200nm. 
This seems to show that the biaxial anisotropy approach may not be enough to solve the 
gate error rate problem and scale MQCA logic to smaller sizes at room temperature. 

APPENDIX A: PROPERTIES OF THE DISCRETIZED STOCHASTIC LLC 
EQUATIONS 

In this appendix we show some of the details of the numerical approach used to solve 
the stochastic LLC equations in the presence of thermal fields. As mentioned before, we 



follow essentially the approach presented in 



2l| . that uses the midpoint rule to discretize the 



stochastic LLC equations. Here we will show that introducing an extra term in the effective 
field that represents the effects of the biaxial anisotropy does not change the two main 
properties of this technique, namely the unconditional preservation of the magnetization 
magnitude and the consistency of the evolution of the free energy. 
The stochastic LLC equations take the form 

dm« = -m(^) X hg^ dt - m^^ x u rfW« + a m^''> x din^'\ (Al) 

where h^jj includes the biaxial term. Applying the midpoint method corresponds to the 
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following replacements: 



dm« ^ (m« 1 - m«) (A2) 



m« ^ (A3) 

hg,(m« t„) hg, + (A4) 

rfW« ^ (W« 1 - W«) (A5) 



that result in the discretized stochastic LLG equations 

' ■""^' + '""' U-(wg'.-w!:') + 



Since every term on the RHS is of the form (^m^*^]^ + m^*)^ x V, it is clear that the RHS van- 
ishes when scalar multiplied by (m^n+-^ + m^^^V while the LHS becomes (|m^*|]^p — |mWp), 



and so we have that 

|m«,r = |m(:)r, (A7) 

which means the midpoint method unconditionally preserves the magnitude of the magne- 
tization. Note that the form of the term added to the effective field does not affect this 



X V. 



2l| have 



property, since the corresponding term on the RHS is still of the form (m^*:^^ + m, 

Another property that the discretized stochastic LLG equations presented in 
is that the change in the discretized free energy is bounded by the work performed by the 
thermal fields on the magnetization for any finite value of the increment At. Their proof of 
this fact relies on the particular form of the effective field, namely that the free energy is 
an at most quadratic polynomial function of the magnetization. Even though when we add 
the biaxial anisotropy term the free energy has a term of degree 4, the result still holds as 
we show below. First we write the free energy g{m) 

- E heit ■ Tj^i^A^ - + my{l- mj)y + m,(l - m^)z), 

19 



where H is the demagnetization tensor, which includes a term corresponding to the uniaxial 
anisotropy, C^"^^^ is the matrix that encodes the dipole-dipole interaction between nanomag- 
nets (and it is symmetric with respect to i and j.) We want to compute Qn+i — Qni where 
Qn = g{mn). Clearly, this will give us an expression in powers of 5m*^*^ = {nifi-i ~ ^n^)- We 
will keep terms up to order ((5m*^*^)^, since from the LLG equations we can see that 5m'^*^ is 
proportional to (w^*|i — W^j, and (w^*]_i — W^*^j is of order At . With that in mind, 
after some algebra we get 

g^^, - 5n ^ E U^'^ • -A/" - E ^^'^ ■ C^''^ ■ ^^'^ - h2t - hlLa^) • 5m« + (A9) 
+- E • M ■ (5m(^) --EE ^^^"^ ■ C^''^ • ^^'^ - 

i i j^i 



where V = l + 2m!^^m^^^ — 3 diag (^(^ri^*^)^, {rrSy^)\, Note that the term multiplying 

5m*^*) in the first sum is exactly — h^yj(m„) (as it should be). Now we go back to the 
discretized LLG equations that have the form 

m« i-m« = -m« , x (hg,(m» , , i„ + ^)At + v (w« , - W») + a (m» , - m«)) 

(AlO) 

where now h^^j also includes the biaxial term. This equation is of the form ra^n+i ~ — 
-m('^ 1 X A, and so if we scalar multiply both sides by A, the RHS vanishes and we get 

n+ 2 



Now we write h^j^ ^m^*|i,i„ + in terms of h^yj (m^*^) (we drop the time, since the 
field does not have an explicit time dependence). After some more algebra, we get 

(m^ii) • Sm^^^t = h'^jf (m») Ai + At E ^m^'^ • C^''^ • ^^^'^ - (A12) 

-Aam» ■ I -^(3M» - 1) - -fA ■ (5m» 

with M'^*) = diag ((mW)^^ (^i'*)n> (^x^)n)- Now, using this in the expression we computed 
for g„^i — gn, and doing even more algebra, we arrive to 

gn+i - = ^ E (w!:ii - W«) • 5m« - E5m(^) • A1«(m„) • <5m», (A13) 
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with A^W(m„) 
we have finally 



At 



a 



1 + m^*^mW'^. Hence, A^*^*^(m„) is positive semidefinite, and so 



gn+i - Qn < z^(ni^*|i 



— m 




(A14) 



At 



which shows that the change in the discretized free energy is always less than the work done 
by the stochastic field during the time interval At. 
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